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Abstract 

We calculate the large-g expansion of the second moment correla- 
tion length at the first order phase transition point of the g-state Potts 
model in two dimensions both in the ordered and disordered phases 
to order 21 in Xj ^Jq. They coincide with each other to the third term 
of the series but differ a little in higher orders. Numerically the ratio 
of the second moment correlation length in the two phases is not far 
from unity in all region of g > 4. The ratio of the second moment 
correlation length to the standard correlation length in the disordered 
phase is far from unity, which suggests that the second largest and 
smaller eigenvalues of the transfer matrix form a continuum spectrum 
not only in the large-g region but also in all the region of (7 > 4. 
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1 Introdunction 



The g-state Potts model in two dimensions exhibits the first order phase 
transition for g > 4. Many quantities^] was solved exactly at the phase tran- 
sition point, among which the correlation length in the disordered phase 
[|, is included. On the other hand, no analytic result had been known for 
the correlation length ^^"^ in the ordered phase at the phase transition point. 
Here the correlation length means the standard one defined from the ratio of 
the largest to the second largest eigenvalues of the transfer matrix and it is 
often called the true correlation length or the exponential correlation length. 
Janke and KapplerJ^ examined the decay rate of the correlation function by 
the Monte Carlo simulation with large sizes of the lattice for q = 10, 15 and 20 
giving a result that was consistent with = ^^'^^ and they made a conjecture 
that this relation would be exact. Igloi and Car Ion |P] evaluated the largest 
and the second largest eigenvalues of the transfer matrix for g = 9 by the 
density matrix renormalization group technique with a result that supports 
this conjecture. Recently the author^] investigated analytically the large-g 
behavior of the eigenvalues of the transfer matrix. He found that from the 
second largest to the L-th largest eigenvalues with L the size of the lattice 
make a continuum spectrum for the thermodynamic limit both in the ordered 
and disordered phases at least in the large-g region and that = ^^'^^ at 
least to order z^^'^ (i.e., in the first 4 terms of the large-g expansion). 

Here we calculate the large-g expansion of the second moment correla- 
tion length ^2nd at the first order phase transition point both in the ordered 
and disordered phases of this model. The second moment correlation length 
also gives important informations on the spectrum of the eigenvalues of the 
transfer matrix. The large-g expansion of the Potts model in two dimensions 
was calculated to order z^^ {z = l/y/q) for the energy cumulants including 
the specific heat by Bhattacharya, Lacaze and Morel |jTO| using a graphical 



method. Tabata and the author O, 12 1 extended the large-g series to or- 



der z"^^ using the finite lattice method]^, |TB|. They also calculated the 
large-g series for the magnetization cumulants including the magnetic suscep- 
tibility to order z^^. These long series enabled them to present the estimate 
for the numerical values of the cumulants which is orders of magnitude more 
precise than that of Monte Carlo simulations. The series also served to con- 
vince us that the correctness of the conjecture by Bhattacharya, Lacaze and 
Morel on the behavior of the divergent quantities in the limit of g ^ 4. It 



2 



is rather straightforward to apply the finite lattice method to the large-g 
expansion of the second moment correlation length, since the algorithm is 
quite parallel with that of the finite lattice method for the low temperature 
expansion of the second moment correlation length for the Ising model in 
three dimensions [1171 . The situation is in contrast to the case of the standard 



correlation length where the continuum spectrum of the eigenvalues of the 
transfer matrix prevented us from applying the method used to obtain the 
low temperature series for the Ising model in three dimensions |ll8|, The 



results of the large-g expansion for the second moment correlation length was 
reported breafiy in reference H and here we describe them in detail. 



2 Expansion series by the finite lattice method 

The model is defined on the x Ly rectangular lattice by the partition 
function 

Z = 5] exp J /5 5] 5s,,s, + + Tia^i + l2y^ + WD^S,,! \ , (1) 

where the spin variable Si at each site i takes the values 1,2, 
represents the pair of nearest neighbor sites and = (xj, Ui) is the coordinate 
of the site i. The phase transition point (3t for /i = 71 = 72 = r/ = is given 
by exp — 1 = y/g. The fixed boundary condition should be taken for the 
ordered phase in which all the spins outside the x Ly lattice are fixed 
to be {si = +1}, and the free boundary condition should be taken for the 
disordered phase. 

The second moment correlation length squared is defined by 

where /i2 is the second moment of the correlation function 

/i2= lim (L^Ly)-^^(ri -rj)2(5,^^i5,^,_i)^, (3) 

fiQ is the zeroth moment of the correlation function, i.e., the magnetic suscep- 
tibility, and d{= 2) is the dimensionality of the lattice. The second moment 
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jj2 can be obtained by the derivative of the free energy density as 



X 2 ( - 

y dhdrj d'fi d-^l 



\nZ{P, /i, 77, 71,72) 



(4) 



/i=»?=7i=72=0 



The algorithm of the finite lattice method to generate the large-g expan- 
sion series for the second moment /i2 is the following. We define H{lx, ly) for 



each Ix X ly lattice (1 < < L^-, ^ <ly < Ly) as 



H{lx, ly) = 2 



92 



92 



92 



dhdrj d'-fi d'-f2 



liaZ{lx,ly) 



(5) 



h=77=7i=72=0 



where Z[lx, ly) is the partition function for the l^ x ly lattice with the fixed and 
free boundary condition for the ordered and disordered phase, respectively, 
and define W of each lattice recursively as 



W{lx,ly) = H{lx,ly) 



i<t'^<lx,l<l'y<ly 



(6) 



X and ly but not on 



Note that H{lx,ly) and W{lx,ly) depend on the size /; 
the position of the origin of the coordinate. Then the second moment of the 
correlation function in the thermodynamic limit is given by 



fi2 = Y. ^(^^^ h 



(7) 



We can prove|]T5[ that the Taylor expansion of W{lx, ly) with respect to 
z = 1/ y/q includes the contribution from all the clusters of polymers in the 



standard cluster expansion that can be embedded into the Ix x ly lattice 
but cannot be embedded into any of its rectangular l'^ x sub-lattices with 
I'x < ^x,ly < ly It is straightforward to understand from the discussion for 



the large-g expansion of the magnetic susceptibility in reference jl^ that the 



I'Xy ly) 



starts from the order of z"' with n = L 



+ ly in 



series expansion of W{1 
the ordered phase and n = Ix + ly — 2 in the disordered phase, respectively. 
So in order to obtain the expansion series to order z^, we should take all the 
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finite-size rectangular lattices for the summation in equation (J^) that satisfy 
Ix + ly < N in the ordered phase and lx + ly — 2<Nm the disordered phase, 
respectively. 

Using these algorithm of the finite lattice method we have calculated the 
series for the second moment H2 at the first order phase transition point to 
order N = 21 in z both in the ordered and disordered phases as 

f^2 = Yl (8) 

n 

The coefficients of the series are listed in table 1. We have checked that all 
of W{lx, ly) with Ix + ly — 2 < N for the ordered phase and with + ly < N 
for the disordered phase have the correct order in z as described above. 

Combining with the large-g series of the magnetic susceptibility given in 
reference [0, we obtain the series for the second moment correlation length 
squared ^1^^ as 

eL = E^n;^"' (9) 

n 

which are listed in table 2. The obtained expansion coefficients for the or- 
dered and disordered phases coincide with each other to order z^ and differ 
a bit from each other in higher orders. 



3 Analysis of the series 



It is known [20 1 that in the limit of the large correlation length 



and , i:r=ic,^(o/eij 



\3 



(10) 



with = -ln(Aj/Ao) and cj =< Ao\0\Aj > where Aj{j = 0, 1,2, ■■■) 
are the eigenvalues of the transfer matrix with Aq the largest one, \Aj > 
is the eigenstate of the transfer matrix corresponding to the eigenvalue A^ 
and O = J2i^si,i with the summation for i running over the sites with a 
fixed coordinate of i/i. The ^1 = — ln(Ai/Ao) is the standard correlation 
length. Later on we will call the eigenstate Ai the first excited state and 
Aj{i = 2, 3, ■ ■ ■) the higher excited states. 
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Table 1: The large-r/ expansion coefficients a„ for the second moment //2 of 
the correlation function aX jS = jSt in the ordered and disordered phases. 



n 


a„ (ordered) 


a„ (disordered) 











1 








2 








3 


4 


4 


4 


44 


44 


5 


344 


336 


6 


2140 


2064 


7 


11676 


11164 


8 


57440 


54652 


9 


261884 


248456 


10 


1121092 


1062320 


11 


4560648 


4320948 


12 


17761800 


16840044 


13 


66657396 


63279396 


14 


242180712 


230307200 


15 


855255360 


815001104 


16 


2945011236 


2812875600 


17 


9914462596 


9493154576 


18 


32704960680 


31397031360 


19 


105910040516 


101949411404 


20 


337252952816 


325538033480 


21 


1057475140292 


1023600956088 



The standard correlation length ^} ' at the phase transition point in the 
disordered phase is known exactly with its asymptotic behavior for — > 4 as 

where 2 cosh 6' = ^ {6 ~ -^q — 4/2 for g — > 4). It is quite natural to expect 
that all of the ^/s behave hke 

^M)^^M)exp(^|!j (12) 
6 



Table 2: The large-g expansion coefficients 6„ for the second- moment cor- 
relation length squared ^2nd at /3 = /3j in the ordered and disordered phases. 



n 6„ (ordered) 6„ (disordered) 












1 


1 


1 


2 


7 


7 


3 


35 


35 


4 


154 


157 


5 


611 


635 


6 


2237 


2363 


7 


7720 


8271 


8 


25461 


27569 


9 


80852 


88271 


10 


249212 


273775 


11 


748239 


825956 


12 


2198991 


2436937 


13 


6335396 


7044039 


14 


17961597 


20027792 


15 


50129172 


56038776 


16 


138121342 


154746684 


17 


375778895 


421876179 


18 


1011869143 


1137902174 


19 


2695052762 


3035669883 



for g — > 4 both in the ordered and disordered phases, in which case 

e£'f ~A(-^)expf^y (13) 



To check the validity of this conjecture we follow the method in reference ||Tl 
adopted to analyze the large-g series of the energy and magnetization cu- 
mulants including the specific heat and the magnetic susceptibility. The 
method was used to convince the validity of the Bhattacharya-Lacaze-Morel 
conjecture on the asymptotic behavior of these quantities for g — > 4. The 
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(a) 



(b) 





Figure 1: The Fade approximants of the second moment correlation length 
squared multiplied by the p-th power of the latent heat (a) in the disordered 
phase and (b) in the ordered phase. 



latent heat L is known exactly with the asymptotic form for g ^ 4 as||2l|| 



TT' 



£~3.exp -- 



(14) 



so, if ^2nd'^ asymptotic form in equation (pTSD, the product ^2nd^^ 

a smooth function of z for g ^ 4 when p = 4, and the Fade approximants of 
the large-g series for this product are expected to converge for p = 4. The 
results are given in figure 1. We can see that the convergence is really good 
around p = 4 both in the ordered and disordered phases. 

The estimated values of the second moment correlation length by setting 
p = 4 are listed in table. 3. We note that the ratio of the estimated values 
of the second moment correlation length in the ordered to disordered phases 
'C2n(i (ordered) /^2nd (disordered) is not far from unity. The ratio is 0.935(5) 
even for g — 4. 

Another interesting and important quantity is the ratio of the second 
moment correlation length C,2nd to the standard correlation length ^i. From 
equation ([T0|) we know that the ratio ^2nd/^i should be less than unity in the 
limit of the large correlation length. If the higher excited states {i = 2,3, ■ ■ ■) 
did not contribute so much, this ratio would be close to unity, as is the case 
in the Ising model on the simple cubic lattice, where ^2nd/C,i = 0.970(5) at 
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Table 3: Estimates of the second moment correlation length ^2nd for each q 
in the ordered and disordered phases. 



Q 


^2nd(ordered) 


^2nd (disordered) 


5 


1261(4) 


1349(6) 


6 


83.1(1) 


88.5(1) 


7 


25.89(1) 


27.45(2) 


8 


13.140(2) 


13.855(5) 


9 


8.3440(5) 


8.751(4) 


10 


5.9965(2) 


6.260(1) 


12 


3.79944(3) 


3.9355(1) 


15 


2.480927(4) 


2.54815(1) 


20 


1.6376632(4) 


1.667680(1) 


30 


1.0633984(1) 


1.0742339(1) 



the critical point in the disordered phase pQ, 22 



In figure 2 we plot the ratio ^2nd/^i for the Potts model in the disordered 
phase. We use the value of ^2nd estimated above from the Pade analysis of 
Cind^^ and the exact value of ^i. The ratio is much less than unity in the re- 
gion of q where the correlation length is large enough. It approaches 0.505(6) 
for g — > 4. This result implies that the contribution of the higher excited 
states is important in the Potts model in two dimensions. It is consistent 
with the fact that the eigenvalue of the transfer matrix for the first excited 
state locates at the edge of the continuum spectrum of the eigenvalues of the 
transfer matrix at least in the large-g region and it strongly suggests that 
these excited states {i = 1, 2, ■ ■ ■ , L^) form a continuum spectrum not only 
in the large-g region but also even in the limit of g — * 4. 



4 Summary 

We calculated the large-g expansion of the second moment correlation length 
in the ordered and disordered phases of the g-state Potts model in two dimen- 
sions and found that they coincide with each other to the third term of the 
expansion but differ a little from each other in higher orders. This suggests 
that, although the first few terms of the large-g expansion for the standard 
correlation length really coincide in the ordered and disordered phases, their 
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0.7 




0.4 I ' ' ' 1 

4 8 12 16 20 

q 

Figure 2: The ratio ^2nd/^i of the second moment correlation length to the 
standard correlation length in the disordered phase. 

higher order terms might be different from each other. Note, however, that 
it is not conclusive, since the second moment correlation length depends not 
only on the spectrum of the eigenvalues of the transfer matrix but also on 
the overlapping amplitude Cj which appeared in equation (|10[). Numerically 
the ratio ^ind/^ind from unity even in the limit of g — 4. We also 

found that ^^nd/^f^ is far from unity for all region of g > 4. It implies that 
higher excited states give significant contributions as well as the first excited 
state and it strongly suggests that these excited states form a continuum 
spectrum (i.e., there is no particle state in the language of the field theory) 
not only in the large-g region but also in all the region of g > 4. Finally we 
point out that it is worthwhile to reanalyze the previous Monte Carlo data 
for the correlation function on the assumption that this would be the case. 
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